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ABSTRACT 

Context. Multiwavelength observations of supernova remnants can be explained within the framework of the diffusive shock accel- 
eration theory, which allows effective conversion of the explosion energy into cosmic rays. Although the models of nonlinear shocks 
describe reasonably well the nonthermal component of emission, certain issues, including the heating of the thermal plasma and the 
related X-ray emission, remain still open. 

Aims. To discuss how the evolution and structure of supernova remnants is affected by strong particle acceleration at the forward 
shock. 

Methods. Analytical estimates combined with detailed discussion of the physical processes. 

Results. The overall dynamics is shown to be relatively insensitive to the amount of particle acceleration, but the post-shock gas tem- 
perature can be reduced by the acceleration to a multiple, even as small as six times, the ambient temperature with a weak dependence 
on the shock speed. This is in marked contrast to models with no particle acceleration where the post-shock temperature is insensitive 
to the ambient temperature and is determined by the square of the shock speed. It thus appears to be possible to suppress effectively 
thermal X-ray emission from remnants by strong particle acceleration. This might provide a clue for understanding the lack of thermal 
X-rays from the TeV bright supernova remnant RX J1713.7-3946. 
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1. Introduction 



2. Dynamics of supernova remnants 



The purpose of this article is to discuss how the evolution of 
supernova remnants may be modified by strong particle accel- 
eration at the forward shock and the resulting implications for 
observational diagnostics, in particular those using X-rays. This 
is an interesting problem in its own right, but is becoming of 
increasing importance because of the rapidly improving obser- 
vational situation; there are now a number of shell-type SNRs 
detected as sources of TeV gamma rays (see for a recent review 
Gabici 2008 and references therein), which certainly proves the 
acceleration of particles to energies of at least 10 14 eV, in partic- 
ular in SNR RX J1713.7-3946 ( |Aharonian efaLpOOTl l, but the 
interpretation of the observations remains controversial, in part 
because of the neglect or misunderstanding of some of these ef- 
fects. In particular, recently a strong claim has been made against 
the hadronic origin of the TeV-gamma-ray emission from the 
young SNRRXJ1713.7-3946 ( |Katz and Waxman |2008[ ), based 
on the lack of thermal X-rays from this TeV-bright object (Slane 
et al. |19991|Cassam-Chenai et al. |2004| |Takahashi et aLT2 008 ; 
Tanaka et al. 2008). However, given the complex nature of the 



problem, which involves the physics of nonlinear shocks as well 
uncertainties related to the heating mechanism of electrons, such 
a claim is premature and further studies are needed. In this pa- 
per we discuss the impact of the efficient particle acceleration 
in strong nonlinear shocks on the heating of the ion component 
of plasma. Our intent is to give a physically motivated account 
of the key processes. We feel that this is more important at this 
stage than the construction of complicated numerical models, 
although these will clearly be required for the detailed under- 
standing of specific sources. 



In essence a supernova remnant, once it has expanded out of the 
initial free-expansion phase where the total energy is dominated 
by the kinetic energy of the expanding ejecta from the explosion 
of the progenitor star, is just a hot bubble of high-pressure gas 
pushing a shell of swept-up material into a low pressure external 
medium. It is important to realise that this is a dynamic situation 
and that even if the system appears to evolve in a quasi-steady 
manner (as in the generalised Sedov solutions) this represents a 
dynamical balance between the conversion of internal pressure 
into kinetic energy of expansion by "PdV" work on the one hand 
and the competing conversion of bulk kinetic energy into random 
particle motion, i.e. pressure, in the shock on the other. 

If the radius of the bubble is R(t) and the explosion energy 
is £sn ~ 10 51 erg then to within factors of order unity the to- 
tal pressure within the bubble has to be 3Es^/4nR 3 with only a 
weak dependence on whether this pressure comes from cosmic 
rays, thermal gas or magnetic fields. More precisely, for rela- 
tivistic accelerated particles and tangled magnetic fields the pres- 
sure is 1/3 the energy density, and for non-relativistic thermal 
particles the pressure is 2/3 the energy density. The average in- 
terior pressure is thus between 1 /3 and 2/3 of 3(£sn -Ek)/4jtR 3 
where is the kinetic energy of the remnant expansion, which 
we can take as some roughly constant fraction of the total en- 
ergy £sn- Equating this internal pressure, again within factors of 
order unity, to the ram pressure of the external medium flowing 
into the shock gives the expansion of the remnant, namely 
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or 

— * R ~ 3/2 { — 
dt \ pi) 



with a difFusive flux from behind the shock into the interior 



1/2 



from which the well-know Sedov self-similar relation, 



Po 

immediately follows. The shock velocity V then scales as 
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It is worth noting that the factors of order unity mentioned 
above all go with the factor Esn/Po an d thus affect the final re- 
sult for the radius as a function of time only to the very weak 
power 1/5. It follows immediately that particle acceleration has 
a relatively small effect on the bulk expansion of the remnant and 
that the basic Sedov scaling for the shock radius as a function 
of time, explosion energy and external density is very robust. 
However the shock structure itself, and the internal structure of 
the remnant, are much more sensitive to the effects of particle 
acceleration. 

The primary effect of putting energy into particle accelera- 
tion is to reduce the amount of energy available for gas heat- 



ing (Ellison 



2000 Blasi et al. 2005), but it also leads to in- 



creased shock compression ratios. In part this is due to the softer 
equation of state associated with the relativistic particles, but the 
more significant reason is the need to constantly pump energy 
into particles near the maximum energy (upper cut-off) of the 
spectrum. There are three distinct reasons for this. 

- (i) One is to compensate for high-energy particles physically 
escaping from the system. This is clearly important if the ef- 
fective magnetic field, and thus the upper cut-off energy for 
the spectrum of accelerated particles, is a decreasing func- 
tion of time as is the case with the currently fashionable 
field amplification models (see e.g., Bell 2004). The effect 
has perhaps been somewhat over-emphasised by the wide- 
spread use of so-called "free escape boundaries" as technical 
devices to regularise steady planar models of shock acceler- 
ation (see, e.g., |Elhson| 1985] |Reville, Kirk and Duffy|2008) . 

- (ii) The second, often overlooked, reason is to compensate 
for the geometrical expansion of the acceleration region as 
the shock evolves. In an expanding and decelerating spheri- 
cal shock the acceleration volume scales as the surface area 
of the shock times the diffusion length-scale, both increasing 
functions of time, and thus an ever larger number of parti- 
cles are needed at a given energy to maintain the particle 
intensity. At energies well below the cut-off the acceleration 
time-scales are very short compared to the dynamical evolu- 
tion times and the system can easily stay in equilibrium, but 
as we approach the upper cut-off energy the time-scales be- 
come comparable and this dilution effect becomes a signifi- 
cant drain on the available supply of accelerated particles. In 
fact, if no other effects (radiative losses, plasma processes, 
collisions etc.) intervene, it is precisely the on-set of this di- 
lution that determines the location of the upper spectral cut- 
off. 

- (iii) Thirdly, as we will see, in the approximately self-similar 
Sedov-like phase of the evolution the interior pressure in the 
remnant is only about half the immediate post-shock value. 
There is thus a gradient in the accelerated particle pressure in 
the down-stream post-shock region which must be associated 



of the remnant (Drury et al. 1989 1. This effect is entirely 



absent in steady planar solutions. 

Compared to a stationary planar shock solution, these effects all 
represent additional energy sinks (from the shock itself; only the 
first is a real energy loss from the whole system) and just as in 
the case of radiatively cooled gas shocks lead to increased shock 
compression ratios (typically 10 or more is seen in simulations). 
This is the fundamental reason why steady one-dimensional 
shock models, such as the original two-fluid models ( Drury and 
Volk 1981)1 and their descendants, are of limited applicability to 



supernova remnants; in addition to the energy loss by advection 
downstream, there has to be a significant flux of energy to par- 
ticles at the upper cut-off O and thus the mass, momentum and 
energy conservation relations for the shock take the approximate 
modified form 



polio = p\U\ = A 
p U 2 = p x U\+P 



(5) 



\p\U\ + U x (8 + P) + 0> 



where Uq is the upstream velocity, U\ the post-shock deceler- 
ated velocity, 8 is the total post-shock energy density of thermal 
and non-thermal particles, P the associated total pressure and we 
further assume that these are negligible upstream of the shock. It 
is important to note that these relations are only approximate for 
an expanding spherical system because the shock is not planar 
and is not in a stationary equilibrium, but the additional energy 
flux term captures much of this non-planar and non-equilibrium 
character. These imply 



A(ffo-tfi) = P 
-P(U + Ui) = Ui (8 + P) + <5 



1 



(6) 
(7) 



and the latter implies that the shock compression ratio is given 
by 



U Q , 28 2<t 

S = = 1 + H . 

Ui P UiP 



(8) 



For non-relativistic particles 8/P = 3/2 rising to 8/P = 3 for 
relativistic particles; thus 1 + 28/ P lies between 4 and 7 with an 
additional positive contribution from the 2^/U\P term. 

We can estimate the order of magnitude of the <t> term as 
follows. From basic shock acceleration theory (see Appendix A) 
the number flux of particles being accelerated upwards in energy 
at momentum p is 



4np : 



f(p)(U -Ui) 



and thus the associated energy flux at the upper cut-off is 

—^flPmadiUo-Ui) CPmax- 



(9) 



(10) 



But this is essentially just the integrand in the total non-thermal 
particle pressure when expressed as an integral over logarithmic 
momentum, 



J ^-4np 2 f(p)dp = j 



\np- 



pvf(p)d\n(p) 



(11) 
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and thus, if the pressure is contributed more or less uniformly 
per logarithmic interval as in the test-particle case, f(p) oc p~ 4 , 
we have that 



Pc(U - uo 



(12) 



of equation of state or particle acceleration. It follows trivially 
that the total interior energy of the remnant (hot gas and cosmic 
rays) is between 3/2 (gas dominated) and 3 (cosmic ray dom- 
inated) times the kinetic energy of the remnant and thus in the 
two extreme cases we have, if gas dominated, 



where A > 1 is the ratio of the total non-thermal pressure to that 
part contributed by particles at the upper cut-off. For an equal en- 
ergy per logarithmic interval spectrum,/(p) oc p~ 4 , it is simply 
the logarithmic range of the spectrum, A = \n(p max /mc) drop- 
ping to A = 2 for Malkov's universal strongly modified shock 
spectrum, f(p) oc p -3 5 , where most of the energy is carried near 
the upper cut-off (Malk ov |1999| l. It is then easy to see that the 
shock compression is given by 



5=1 + 



2<5 



P-2P C /A' 



(13) 



which can easily be quite large for significant particle accelera- 
tion and small values of A. Writing Pc for the thermal pressure 
and P c for the non-thermal pressure and assuming the latter to 
be dominated by relativistic particles we have 



1 + 



3P G + 6P C 
Pg + PcQ-W 



(14) 



It is very remarkable that even at this crude macroscopic level 
we can see a connection between extreme acceleration efficiency 
with the shock compression tending to infinity and a spectral 
hardening towards an exponent of 3.5 corresponding to A = 2 as 
required by Malkov's asymptotic solution. 

This large compression ratio allows us to retrospectively jus- 
tify the simplified relations we started with using an approxima- 
tion originally due to |Chernyi (1957]) (see also Zel'dovich and 
Rais er |1966[ ). The key is to assume that all the swept-up matter 
is concentrated in a thin shell immediately behind the shock with 
radial velocity Uq - U\ and total mass 



M = yflV 



(15) 



If the pressure in the interior of the remnant is P m , then 
Newton's law of motion applied to an element of the shell gives 



dt 



[M(U - Ui)] = 4nR 2 P inl . 



(16) 



For a self-similar solution the interior pressure Pj nt will be some 
fixed fraction a of the post-shock pressure A(Uq - U{). Thus 



d_ 
dt 



[M(U - Ui)] = 4ttR 2 Pi) U (U - U x )a 



or, noting that Uq = dR/dt, 



R d(Uo-Ui) 



u -u. 



dR 



3(a- 1) 



(17) 



(18) 



from which it is easy to see that the self-similar Sedov-like so- 
lution with R oc t 2/s and U oc R^l 2 requires a — 1/2. Thus the 
interior pressure is half the immediate post-shock pressure. That 
there has to be such a pressure gradient in the shell is obvious 
because the material that was shocked at early times is mov- 
ing too fast and would overtake the shell were it not decelerated 
by an adverse pressure gradient directed towards the interior of 
the remnant. It is however remarkable that the Sedov scaling re- 
quires the total interior pressure to be half the ram pressure of 
the shock, again with only a very weak dependence on details 



Pint 



-SN> 



5 

2 3£sn 
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and if cosmic ray dominated, 



Ek 
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(21) 
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Strictly speaking, in the gas dominated case we should not ne- 
glect the thickness of the shell, and we should also allow for the 
fact that the post-shock gas velocity is only 3/4 of the shock 
velocity. Also in the cosmic-ray dominated case, if O has a sig- 
nificant component from genuine escape as distinct from geo- 
metrical dilution, this will slightly affect the Sedov exponent. 
However the main point of this section is to demonstrate that the 
interior pressure is rather tightly constrained and is of order 0.4 
to 0.25 times the explosion energy divided by the remnant vol- 
ume, and that the shock expansion is well approximated by the 
standard Sedov formula. 

Because the total pressure is essentially fixed, if we imagine 
more and more efficient particle acceleration putting more pres- 
sure into Pc, the thermal gas pressure Pc has to decrease. At 
the same time the nonlinear reaction terms are making the shock 
more compressive and increasing the downstream density. Both 
effect cause the post-shock temperature Ti oc P G /p to decrease so 
that we expect substantially lower post-shock gas temperatures 
than in pure gas models of SNRs. Once a fluid element (in this 
case a spherical annulus) has been shocked it gradually moves 
back through the dense shell expanding and dropping in pres- 
sure until it merges into the tenuous interior region. The shell 
thickness is determined by 



4nR 2 ARp 



4n 



-R 3 p 



(23) 



or AR « R/3s 
shell is thus 



R/12. The time taken to transit through the 



At 



AR 



R r 
3r~k 



1 R 

3k 



(24) 



for a Sedov scaling, close enough to t. 

Thus if we consider a fluid element that is shocked at time 
to, it exits the shell at approximately (1 + 5/6)fo with a pressure 
half the post-shock pressure at that time, which in turn is a fac- 
tor (1 + 5/6)~ 6/5 ~ 1/2 what it was at the initial time to- Thus 
between fo and about 2to the pressure has to drop by a factor of 
about 4. After the fluid element has entered the interior region 
it is a reasonably good approximation to say that it then simply 
expands in pressure equilibrium with the interior of the remnant 
where the total pressure is dropping as 



Pin oc E SN R 3 cc r 6/5 . 



(25) 



Physically this is because the interior flow is subsonic and any 
pressure variations get equilibrated on the sound-crossing time 
which is short compared to the dynamical time of the remnant. 
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So the history of the parcel of shocked gas is of an initial, rather 
rapid, drop in pressure by a factor of roughly four between fo and 
2to followed by a slower power-law r 6 ^ 5 decline as the remnant 
expands. 

Neglecting for the moment any diffusion of the accelerated 
particles out of the shell and splitting the total pressure into a 
"gas" component Pc from thermal ions and a "cosmic ray" com- 
ponent Pc from the relativistic accelerated particles, we have 

(26) 
(27) 



The pressure per logarithmic interval of momentum is given 
by the integrand in 



P 



5/3 



and thus the gas pressure drops faster than the cosmic ray pres- 
sure in the expansion. It follows that if the two were initially 
comparable, the cosmic ray pressure will dominate at later stages 
of the expansion and the density will have to drop according to 

(28) 
(29) 
(30) 



P oc r 9/1 ° 

p c cc r 6 / 5 

-3/2 



(31) 



P C oc t 

and thus the thermal ion temperature will drop as 

rp P G -3/5 
7j oc oc t ', 

P 

significantly faster than in a remnant not dominated by cosmic 
rays where the ion temperature drops as r 12 ^ 25 . This of course 
ignores the diffusion of cosmic rays in the interior, but as long 
as the pressure is more or less uniform this will not affect the 
result significantly. The key point is that the cosmic rays come 
to dominate the pressure as the system expands and this leads to 
a more rapid cooling of the thermal ions. 

In summary therefore the effect of increasing particle accel- 
eration is to leave the rate of expansion of the remnant, and 
the Sedov scaling with external density and explosion energy, 
largely unchanged but to make the shock more compressive, the 
thermal gas colder and the internal pressure become cosmic -ray 
dominated. 



2.1. Self-regulated acceleration 

The accelerated particle pressure is determined by the interaction 
of two processes: 

- (i) The first is the injection of particles into the acceleration 
process at low energies (the injection momentum /?j n jfor pro- 
tons is typically taken as a few times the thermal energy of 
the shock-heated protons) which determines the number of 
particles available for acceleration. There are good reasons 
for thinking that the injection preferentially favours high 
rigidity species, but we have no soundly based theory which 
would allow us to quantitatively predict the injection of vari- 
ous species. However at least for non-relativistic shocks the- 
oretical considerations and the evidence from heliospheric 
shocks points to a relatively easy injection of ions as long as 
the shock is of reasonable strength. 

- (ii) The second is the separation between the injection scale 
and the upper cut-off of the accelerated particle spectrum. 
The key difference between supernova remnant shocks and 
heliospheric shocks is that in the former the length and time 
scales are such that a shock-accelerated spectrum can be es- 
tablished which stretches over many decades giving a very 
large lever-arm between injection and the upper cut-off mo- 
mentum p mdx . For the power-law spectra predicted by the lin- 
ear theory of shock acceleration this large lever arm makes 
the total pressure very sensitive to the spectral slope. 



r 

•J Pin 



— Anp f(p)p— 
3 P 



(32) 



which scales as p f(p) in the non-relativistic region and as 
p 4 f(p) in the relativistic region of the spectrum. Thus if the ac- 
celerated particle spectrum is harder than p~ 5 the pressure inte- 
grand is a rising function of momentum in the non-relativistic 
region, and if it is harder than p~ A it is a rising function in the 
relativistic region. As is well know simple test particle theory 
gives a slope for the power-law spectrum of accelerated parti- 
cles which depends only on the compression ratio of the shock, 



pdf 
fdp 



3 s 



1 



(33) 



Thus, for a strong shock in a ideal monatomic gas (i.e. an adi- 
abatic exponent of 5/3) with compression ratio s = 4 we get a 
power-law index of 4, so that the pressure rises linearly with mo- 
mentum in the non-relativistic region and is then logarithmically 
divergent in the relativistic region. 

The effect of the reaction of the accelerated particles back on 
the flow is twofold (see Malkov and Drury 2001 for a review). 



- (i) It makes the overall shock structure more compressive, 
and thus for particles that are sufficiently energetic to see 
this increased shock compression the spectrum hardens with 
the consequence that the pressure integrand has a maximum 
close to the upper cut-off. 

- (ii) At low energies however the subshock, where the dissi- 
pation is determined by collisionless shock physics operating 
on plasma scales, is weakened and the associated spectrum 
of accelerated particles becomes softer. 

Thus the overall spectrum becomes concave, with a softer por- 
tion at low energies becoming harder at high energies and then 
cutting-off, in an approximately exponential fashion, at some 
maximum energy. 

A number of approaches have been used to calculate such 
modified spectra involving various approximations, and the re- 
sults are in good general agreement (although it must be said that 
all these approaches make the questionable assumption, explic- 
itly or implicitly, that the shock structure on intermediate scales 
is quasi-stationary). The fundamental model behind these solu- 
tions is that the shock modification has to adjust itself to throt- 
tle back the injection and initial acceleration to avoid an exces- 
sive energy demand, and this requires the modification to signif- 
icantly weaken the subshock compression. In the limit that that 
injection energy can be treated as small compared to the proton 
rest mass energy, so that there is a long lever-arm between injec- 
tion and the relativistic regime, the above arguments suggest that 
the spectrum in the non-relativistic region should be close to p~ 5 
and thus the subshock compression should approach the critical 
value of i su b ~ 2.5. 

If the subshock compression is weaker than 2.5 the pressure 
integrand can have no local minimum separating a "thermal" 
population from a "non-thermal" component and instead there 
is just a very weak non-thermal tail on a dominant thermal pop- 
ulation. On the other hand, if the subshock is stronger than com- 
pression 2.5, the pressure integrand will start to rise, there will 
be a minimum in the pressure integrand, and we can usefully 
talk about two particle populations; a bulk "thermal" component 
and a non-thermal component where the energy is concentrated 
in relativistic particles. 
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This simple picture is confirmed by the nonlinear models, see 
e.g. the discussion in Berezhko and Ellison ( 1999) and in partic- 
ular their figure 6(b), which show that for efficient injection and 
strong shocks the subshock compression is typically between 2 
and 3 with a rather weak dependence on total Mach number. 
The recent study ( [Kang et al. ||20 08) also finds an asymptotic 
subshock compression of 3.2(M/10) 4 for M > 10 where M is 
the shock Mach number. 



3. Ion temperatures 

If the subshock compression automatically adjusts to a value 
close to 2.5 to 3, it follows that the sub-shock Mach number, 
and the amount of gas heating in the subshock, is essentially 
fixed. This has the very important consequence that the down- 
stream gas temperature is determined largely by the upstream 
gas temperature and not, as in standard shock models, by the 
shock speed. The inflowing gas is first heated by adiabatic com- 
pression in the shock precursor from the far upstream value To 
to a value T\ just in front of the subshock of 



(34) 



where s pie is the precursor compression and then on passage 
through the compression 2.5 subshock by an amount 

T 2 = y^i (35) 

(assuming a y = 5/3 polytropic equation of state for the gas). 
This rises to 



11 



(36) 



for a subshock of compression ratio 3, still a relatively modest 
increase (if there is significant wave dissipation in the precursor 
it is of course possible to raise the gas temperature to higher 
levels ( |Volketal. \\9M\ ). 

If the total compression, as suggested by numerous simula- 
tion studies, is of order 10, with a factor of 4 in the precursor and 
2.5 in the subshock it follows that the downstream temperature 
is just the far upstream temperature multiplied by a factor of 



4 2/3 



12 



6.05 



(37) 



and is independent of the shock speed. This is in marked con- 
trast to the standard picture where the downstream temperature 
is determined purely by the shock speed (with a strong quadratic 
dependence) and is only very weakly dependent on the upstream 
temperature according to, 



:kT G 



1 



-U 



(38) 



k being Boltzmann's constant and m p the proton mass. In reality 
the total compression does increase with shock Mach number 
( |Kang et al. [2008 suggest roughly as the one third power, al- 
though this calculation includes Alfven wave heating effects; a 
naive estimate without wave heating suggests that the exponent 
might be as high as 3/4) and thus there is a weak dependence of 
the post-shock temperature on shock speed (scaling somewhere 
between 2/9 and 1/2 power) but nothing as strong as in the un- 
modified case. 

Clearly this is an extreme limit, but the rather surprising an- 
swer to the question, how cold can the postshock gas be in the 



limit of strong particle acceleration, is as cold as six times the 
upstream temperature! 

Motivated by the above discussion we now give a simple de- 
scription of the gas compression and heating in a supernova rem- 
nant with strong particle acceleration. As an illustrative example 
we follow a fluid element as it is swept up by the shock and 
compressed, initially to a density of 10 times ambient and tem- 
perature of order 6 times ambient at remnant age to. We assume 
that it then expands and cools so that the total pressure drops ini- 
tially by a factor of order 4 over the time to < t < 2to and then 
in pressure equilibrium with the interior. Consistent with the as- 
sumption of strong acceleration we assume that the cosmic ray 
pressure dominates the thermal gas pressure (note that if this is 
true initially, it is true at all later times). In this situation we have 
P c oc p 4 / 3 and T c cc P G /p cc p 2 / 3 so that T c cc p\l 2 . 

Thus the gas temperature can be approximated as 



Tc 
To 



6 — 3| 1 

\to 



-3/5 

and the density as 



to < t < 2/ 
t > 2to 



6T 



Po 



3/2 



(39) 
(40) 

(41) 



This is illustrated in Fig.[T] 



Temperature history 
T 



E? 




Fig. 1. The approximate temperature history as a function of 
time t of an element of gas shocked at time to- The graph shows 
a sketch of the plasma temperature, T c , in terms of the far up- 
stream ambient temperature To for case of no wave heating in the 
precursor, strong cosmic ray acceleration and an assumed total 
compression of 10. 



Even assuming full thermal equilibrium between the elec- 
trons and the ions, it is clear that the gas temperature can be 
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reduced to the point where thermal X-ray emission is severely 
suppressed unless the upstream medium is itself already very 
hot. It is interesting in this context to note the observational sug- 
gestion for substantially lower than expected ion and electron 
temperatures seen in the Magellanic cloud remnant IE 0102.2- 
7219 ( [Hughes et al."|2000) . 



4. Conclusions 

We have attempted to outline the physics behind the modifica- 
tions of SNR structure and dynamics by particle acceleration. 
While the overall dynamics are relatively insensitive, we have 
shown that the gas heating can in principle be strongly sup- 
pressed to the point where thermal X-ray emission is no longer 
expected. This has important consequences for the on-going de- 
bate on the origin of the TeV emission detected from a number 
of young shell-type SNRs. 
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Appendix A: The accelerated particle flux 

The expression (9) for the number flux of particles being accel- 
erated upwards in momentum at the shock can be derived in a 
number of different ways and is fundamental to the theory of 
shock acceleration. 

Perhaps one of the more physical derivations which empha- 
sises the generality of the result is the following. We begin by 
considering a general oblique non-relativistic MHD shock in the 
standard de Hoffman-Teller frame where the plasma velocity is 
always parallel to the magnetic field and the electric field van- 
ishes identically. Let the plasma velocity be U i upstream and U 2 
downstream with both velocities small relative to the speed of 
light c and let the shock normal direction be n. 

Let us now consider a relativistic particle of momentum p 
and velocity v sufficiently energetic that its gyroradius is much 
larger than the shock thickness. In the de Hoffman-Teller frame 
it simply crosses the shock front ballistically with no change 
to momentum or velocity. However we conventionally mea- 
sure particle energies relative to a frame co-moving with the 
local plasma, and in these reference frames there is a change. 



Transforming into the local plasma frame from the de Hoffman- 
Teller frame induces a first-order change in the magnitude of the 
particles momentum of 



Ap 



p U 



(A.l) 



as can be trivially verified by carrying out the Lorentz transfor- 
mation of the particle's energy momentum 4-vector and neglect- 
ing terms of second order and higher in U/c. There is thus a 
change in momentum, measured in the local co-moving plasma, 
of 



Ap = 



p-(Ur- U 2 ) 



(A.2) 



for the particle crossing the shock. 

Let us now consider a distribution of particles incident with 
different velocity vector directions on the shock. If we consider 
velocity vectors within a certain solid angle d£2 the flux across 
the shock front is 



p 2 f(p)v-ndQ 



(A3) 



and thus the flux of particles upwards through momentum level 
p at the shock is 



u 2 ) 



dn. 



(A.4) 



In general this is all that one can say, but if the distribution func- 
tion is close to isotropic there is a further significant simplifi- 
cation. In the special case of isotropic distributions, by symme- 
try the integral must reduce to a scalar times the dot product 
n ■ (Ui — I/2) an d by evaluating in the special case of a parallel 
shock (with 9 the angle between the particle velocity vector and 
the shock normal) we easily find 



¥(P) = P 3 /(P) J n-(U l -U 2 )cos 2 (0)d^ 
-p 3 f(p)n ■ (t/j - U 2 ) 



4n 

y 



(A.5) 



as in equation (9). 

This emphasises that the result depends only on the distribu- 
tion function of the accelerated particles being close to isotropy 
at the shock and the velocity shifts into the local plasma frame 
being sub-relativistic; the two assumptions are of course closely 
related the main problem in the theory of relativistic shock ac- 
celeration being that the distributions are highly anisotropic at 
the shock. 

Another, perhaps more conventional, route to this result is to 
write the basic diffusive transport equation in conservation form. 
Normally the transport equation is written in the form 



|f + t /.v/=iv-i/^ + vuv/), 

at 3 op 

but we can equivalently write this as 

(4np 2 f(p)) + V-{Anp 2 f(p)U-Anp 1 KVf) 



(A.6) 



at 



+ 



dp 



4np- 







(A.7) 



which can be interpreted as stating that the time rate of change 
of the number density, plus the physical space divergence of the 
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advective and diffusive fluxes, plus the momentum space diver- 
gence of the acceleration flux is zero. Thus the acceleration flux 
in momentum space is 

V=^y-f(p)(-V-U) (A.8) 

and is directly proportional to the flow compression. In the case 
of a shock the compression is concentrated in a Dirac delta dis- 
tribution at the shock itself 

-V -U = n-(U 1 -U 2 )S(x) (A.9) 
and thus we get a localised acceleration flux at the shock, 

* = ^f-f(P)n ■ (Ui ~ U 2 )S(x) (A. 10) 

in agreement with equation (9). 



